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By calculation and analysis of the bare conduction bands in a large number of hole-doped 
high-temperature superconductors, we have identified the energy of the so-called axial-orbital as 
the essential, material-dependent parameter. It is uniquely related to the range of the intra-layer 
hopping. It controls the Cu 4s-character, influences the perpendicular hopping, and correlates 
with the observed T c at optimal doping. We explain its dependence on chemical composition and 
structure, and present a generic tight-binding model. 



PACS numbers: 74.25. Jb, 74.62. Bf, 74.62.Fj, 74.72-h 
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The mechanism of high-temperature superconductiv- 
ity (HTSC) in the hole-doped cuprates remains a puz- 
zle Q. Many families with Cu02-layers have been syn- 
thesized and all exhibit a phase diagram with T c going 
through a maximum as a function of doping. The pre- 
vailing explanation is that at low doping, superconduc- 
tivity is destroyed with rising temperature by the loss 
of phase coherence, and at high doping by pair-breaking 
[^|. For the moienate-dependence of T c at optimal dop- 
ing, T cmax , the only known, but not understood, sys- 
tematics is that for materials with multiple Cu02-layers, 
such as HgBa2Ca„_iCu„02n+2, T cmax increases with the 
number of layers, n, until n —3. There is little clue as 
to why for n fixed, T cmax depends strongly on the fam- 
ily, e.g. why for n=l, T cmax is 40 K for La2CuC>4 and 
85 K for T^E^CuOg, although the Neel temperatures 
are fairly similar. A wealth of structural data has been 
obtained, and correlations between structure and T c have 
often been looked for as functions of doping, pressure, 
uniaxial strain, and family. However, the large num- 
ber of structural and compositional parameters makes 
it difficult to find what besides doping controls the su- 
perconductivity. Insight was recently provided by Seo et 
al. H who grew ultrathin epitaxial Lai.gSro.iCuC^ films 
with varying degrees of strain and measured all relevant 
structural parameters and physical properties. For this 
single-layer material it was concluded that the distance 
between the charge reservoir and the Cu02-plane is the 
key structural parameter determining the normal state 
and superconducting properties. 

Most theories of HTSC are based on a Hubbard model 
with one Cuti^.^-like orbital per CuC-2 unit. The one- 
electron part of this model is, in the k-representation: 

e (k) = — 2t (cos k x + cos k y ) + At' cos k x cos k y 

-2t"(coa2k x + cos2k y ) + ..., (1) 

with t,t',t",... denoting the hopping integrals (> 0) 
on the square lattice (Fig.Q). First, only t was taken 
into account, but the consistent results of local-density 
approximation (LDA) band-structure calculations Q 
and angle-resolved photoemission spectroscopy (for over- 



doped, stripe- free materials) have lead to the current 
usage of including also t' , with t' /t —0.1 for La2Cu04 and 
t'/t -0.3 for YBa 2 Cu 3 7 and Bi 2 Sr 2 CaCu 2 8 , whereby 
the constant-energy contours of expression (|l|) become 
rounded squares oriented in respectively the [11]- and 
[10]-directions. It is conceivable that the materials- 
dependence enters the Hamiltonian primarily via its one- 
electron part (0), and that this dependence is captured 
by LDA calculations. But it needs to be filtered out: 
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The LDA band structure of the best known, and only 
stoichiometric optimally doped HTSC, YBa2Cu30y, is 
more complicated than what can be described with the 
t-t' model. Nevertheless, careful analysis has shown £| 
that the Zow-energy, Zayer-relatcd features, which are the 
only generic ones, can be described by a nearest-neighbor, 
tight-binding model with four orbitals per layer (Fig. [l]), 
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hopping tg S proceeding via the diffuse Cu s-orbital whose 
energy e s is several eV above the conduction band. Also 
the intralayer hoppings t' , t", ... beyond nearest neighbors 
in (|l]) proceed via Cus. The constant-energy contours, 
Si (k)=E, of this model could be expressed simply as H: 



1 - u -d(e) + (1 + u)p(e) = 



He) 



(2) 



in terms of the coordinates u = | (cosfc x + cosk y ) and 
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cos k y ) , and the quadratic functions 
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which describe the coupling of O a /i,p x / y to respectively 
Cud x 2_ y 2 and Cus. The term proportional to p(e) in 
(0) describes the admixture of O a /bP z orbitals for dim- 
pled layers and actually extends the four-orbital model 
to a six-orbital one [Q. For e near the middle of the 
conduction band, d (e) , s (e) , and p (e) are positive, and 
the energy dependence of d(e) may be linearized (d>0), 
while that of s (e) and of p (e) may be neglected. p=0 



for flat layers and p—s j (1 



for layers dimpled so 



as to yield extended saddlepoints. The bilayer bonding 
and antibonding subbands have e s - values split by 
Now, if s s were infinitely far above the conduction band, 
or t sp vanishingly small, the right-hand side of ^ would 
vanish, with the result that the constant-energy contours 
would depend only on u. The dispersion of the conduc- 
tion band near the Fermi level would thus be that of the 
one-orbital model (@) witn M 1 ~ P) l^ d ' and t'=t"=0. 
For realistic values of e s and t sp , the conduction band 
attains Cu s-character proportional to v 2 , thus vanishing 
along the nodal direction, k x —k y , and peaking at (7r, 0) , 
where it is of order 10 per cent. The repulsion from the 
Cus-band lowers the energy of the van Hove singulari- 
ties and turns the constant-energy contours towards [10]. 
This same ^-dependence pertains to the interlayer split- 
ting caused by tjr s in a multilayer material. In order to 
go from © to (|), 

1 2r ... 1/2 



1 — u + s 1 — 2ru ' 



with 



1 + s- 



(3) 



was expanded in powers of 2ru. This provided ex- 
plicit expressions, such as: t = [1 — p + o (r)] /4c?, t' — 
[r + o (r)] /4d, and t" = ht' + o (r) , for the hopping inte- 
grals of the one-orbital model in terms of the parameters 
of the four(six)-orbital model and the expansion energy 
~ £f- Note that all intralayer hoppings beyond nearest 
neighbors are expressed in terms of the range-parameter 
r. Although one may think of r as t' /t, this holds only for 
fiat layers and when r<0.2. When r>0.2, the series (|l|) 
must be carried beyond t" . Dimpling is seen not to influ- 
ence the range of the intralayer hopping, but to reduce t 
through admixture of O a /i,p z . In addition, it reduces t v &. 

Here, we shall generalize this analysis to all known fam- 
ilies of HTSC materials using a new muffin-tin-orbital 
(MTO) method [§ which allows us to construct mini- 
mal basis sets for the low-energy part of an LDA band 
structure with sufficient accuracy that we can extract 
the materials dependence. This dependence, we find to 
be contained solely in e s , which is now the energy of the 
axial orbital, a hybrid between Cus, Cud 3i: 2_ 1 , apical- 
oxygen O c p z , and farther orbitals on e.g. La or Hg. The 
range, r, of the intralayer hopping is thus controlled by 
the structure and chemical composition perpendicular to 
the Cu02-layers. It turns out that the materials with the 
larger r (lower £ s ) tend to be those with the higher ob- 
served values of T cmax . In the materials with the highest 
Tcmax, the axial orbital is almost pure Cu4s. 
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FIG. 2. LDA bands calculated with the NMTO method 
in the body-centered tetragonal structure. The dashed band 
was obtained using the Bloch sum of MTOs with N=0 and Cu 
d x 2_ y 2 symmetry at the central site. F (0,0,0) , D (n, 0, 0) , 
Z (2tt, 0, 0) = (0, 0, 2%/c) , X (tt, tt, 0) . 



It should be noted that r describes the shape of the 
non-interacting band in a 1 eV-range around the Fermi 
level, whose accurate position is unknown because we 
make no assumptions about the remaining terms of the 
Hamiltonian, inhomogeneities, stripes, a.s.o. 

Fig. ^ shows the LDA bands for the single-layer ma- 
terials La2Cu04 and Tl 2 Ba2Cu06. Whereas the high- 
energy band structures are complicated and very differ- 
ent, the low-energy conduction bands shown by dashed 
lines contain the generic features. Most notably, the dis- 
persion along TDZ is suppressed for T^L^CuOe rela- 
tively to La2Cu04, whereas the dispersion along TXZ is 
the same. This is the f 2 -effect. The low-energy bands 
were calculated variationally with a single Bloch sum of 
Cu d x 2 _ y 2-like orbitals, constructed to be correct at an 
energy near half-filling. Hence, these bands agree with 
the full band structures to linear order and head to- 
wards the pure Cu d a .2_ J( 2-levels at T and Z, extrapolating 
across a multitude of other bands. This was explained in 
Ref. |@]. Now, the hopping integrals t, t', t" ', .... may be 
obtained by expanding the low-energy band as a Fourier 
series. This yields: £=0.43 eV in both cases, i'/t=0.17 
for La2Cu04 and 0.33 for Tl2Ba2Cu06, and many fur- 
ther inter- and intralayer hopping integrals Q| . 

That all these hopping integrals and their materials- 
dependence can be described with a generalized four- 
orbital model, is conceivable from the appearance of the 
conduction-band orbital for La2Cu04 in the rEz-plane 
(Fig. ||). Starting from the central Cu atom and going in 
the iE-direction, we see 3d x 2_ y 2 antibond to neighboring 
O a 2p x , which itself bonds to 4s and antibonds to 3d 3z 2_ 1 
on the next Cu. From here, and in the z-direction, we see 
4s and 3d 3z 2_ 1 antibond to O c 2p Zl which itself bonds to 
La orbitals, mostly 5d 3z 2_ 1 . In the y-direction, 4s anti- 
bonds and 3d 3z 2_i bonds to Ob 2p y . For Tl2Ba2Cu06, we 
find about the same amount of Cu 3d x 2_ y 2 and O a /b 2p x / y 
character, but more Cu4s, negligible Cu3c? 3z 2_ 1 , much 
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less O c 2p z , and T16s instead of La5d3 2 2_ 1 character. 
That is, in Tl 2 Ba 2 CuOg the axial part is mainly Cu4s. 




FIG. 3. Left: N=0 MTO describing the Cu d^^-like 
conduction band in La2Cu04. The plane is perpendicular to 
the layers and passes through Cu, O a , O c , and La. Right: 
Schematic diagram giving the energy e s of the axial orbital in 
terms of the energies of its constituents and their couplings. 



Calculations with larger basis sets than one MTO per 
Cu02 now confirm that, in order to localize the orbitals 
so much that only nearest-neighbor hoppings are essen- 
tial, one needs to add one orbital, Cu axial, to the three 
standard ones |?J . The corresponding four-orbital Hamil- 
tonian is therefore the one described above in Fig. |l] and 
Eqs. (|^)-(^). Note, that we continue to call the en- 
ergy of the axial orbital e s , and its hopping integral with 
Oa/bPx/y tsp- Calculations with this basis set for many 
different materials show that, of all the parameters, only 
e s varies significantly [Q. This variation can be under- 
stood in terms of the couplings between the constituents 
of the axial orbital sketched in the right-hand panel of 
Fig. H We first form the appropriate O c p z -like 5-atom 
hybrid Cud 3z 2_ 1 - 20 c p z -2 La with the energy (?]] 
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c La 
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(4) 



and then couple this to the Cus-orbital to yield the 
energy e s — e s + 2t 2 sc / [ep — e c ) of the axial orbital. 
Here, the energies of the pure Cus- and O c p z -orbitals 
are denoted e s and e E , respectively, while their hopping 
integral is t sc . The energy of the Cu d 3z 2_ 1 -orbital is 
e z 2, and its hopping integrals to O a / p x / y and O c p z 
are respectively t pz 2 and t cz i. In deriving Eqs. 
we have exploited that t 2 z2 /t 2 p <§; 6 e F ~y and that 

tlJt% « ZZ^X^/l - Although specific for La 2 Cu0 4 , 
Eq. (El) is easy to generalize. 
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FIG. 4. Calculated range parameter, r, for single-layer 
materials vs. the distance (in A) between Cu and O c . The 
lines result from rigid displacements of O c . 

In Fig. ^ we plot the r-values for single-layer mate- 
rials against the distance dcu-o c between Cu and api- 
cal oxygen, r increases with <icu-o c because e s is low- 
ered towards when the coupling between O c p z and 
Cuc?3 2 2_ 1 /s is weakened. Since t cz i oc d^_ and 
t sc oc oIq 2 _ , increasing the distance suppresses the 
Cud 3z 2_ 1 content, which is then important in La2Cu04, 
but negligible in Tl 2 Ba2Cu06 and HgBa 2 CuC>4. This is 
also reflected in the slopes of the lines in Fig. || which 
give r vs. dcu-o c for each material. The strong slope 
for La 2 Cu04 explains the findings of Seo et al. ||, pro- 
vided that r correlates with superconductivity. That the 
Bi-point does not fall on the La-line is an effect of Bi be- 
ing different from La: Bi 6p z couples stronger to O c 2p z 
than does La5d 3z 2_ 1 . The figure shows that upon reach- 
ing HgBa 2 CuC>4, r is saturated, e s ~ e s , and the axial 
orbital is almost purely Cu4s. 
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FIG. 5. Correlation between calculated r and observed 
Icmax. Filled squares: Single-layer materials and most bond- 
ing subband for multilayers. Empty squares: Most antibond- 
ing subband. Half-filled squares: Non-bonding subband. Dot- 
ted lines connect subband-values. Bars give fc z -dispersion of 
r in primitive tetragonal materials, a-m |8H20| 

Fig. Bj hints that for single-layer materials r might cor- 
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relate with the observed T cmax . But the experimental 
uncertainties of both T cmax and the structural parame- 
ters are such that we need better statistics. Therefore, we 
plot the observed T cmax against the calculated r- values 
for nearly all known hole-doped HTSCs in Fig. ||. For 
the single-layer materials, we observe a strong correla- 
tion between r and T cmax , which seems to be contin- 
ued in the bonding subband for the multilayer materi- 
als (filled squares). This indicates that the electrons are 
delocalized over the multilayer [pl|, and that T cmax in- 
creases with the number of layers for the same reason 
that it increases among single-layer materials; the mul- 
tilayer is simply a means of lowering e s further, through 
the formation of Cu4s-Cu4s bonding states. This is con- 
sistent with the celebrated pressure-enhancement 0| of 
T c in HgBa2Ca2Cu308- One might attempt to increase 
7c max , say for YBa2Cu30y, by substituting Y with a 
smaller cation, e.g. Sc. This has not been done, but 



a larger cation, La, was recently inserted |12|, and that 
caused T cmax to drop from 92 K to 50 K. Using the ob- 
served structure of LaBa2Cu307, we have calculated the 
r-values and included them in Fig. ^. Here again, the 
bonding subband is seen to follow the trend! That T cmax 
eventually drops for an increasing number of layers, is 
presumably caused by loss of phase coherence. 

Interlayer coupling in bet La2Cu04 mainly proceeds 
by hopping from O c p z at (0,0, zc) to its four near- 
est neighbors at (±|, ±5, (5 — z) c) , and is therefore 
taken into account by adding to e E on the right-hand 
side of (Q) the term — 8t^r c cos \k x cos ^k y cos \ck z . In 
primitive tetragonal materials, the corresponding term is 
merely oc coscA^ because the Cu02-layers are stacked 
on top of each other, e.g. in HgBa2Cu04, the in- 
terlayer coupling proceeds from O c p z at (0,0, zc) via 
Hg6s/6p 2 at (0,0, c/2) to O c p z at (0, 0, (1 - z) c) . Pe- 
riodic interlayer coupling thus makes e s depend on k z , 
and this passes onto the conduction band a A: z -dispersion 
oc v 2 cos \k x cos \ k y cos \ck z in bet and oc v 2 cos ck z in 
tetragonal structures. Fig. [| shows how the fc 2 -dispersion 
of r decreases with contraction of the axial orbital. 

Our identification of an electronic parameter, r or e s , 
which correlates with the observed T cmax for all known 
types of hole-doped HTSC materials should be a use- 
ful guide for materials synthesis and a key to under- 
standing the mechanism of HTSC. With current k-space 
renormalization-group methods one could for instance 
investigate the effect of the band shape on the leading 
correlation-driven instabilities |23). Moreover, the possi- 
bility that a longer hopping-range leads to better screen- 
ing of the Coulomb repulsion, maybe even to overscreen- 
ing, could be studied. Increased diagonal hopping, t' , 
might lead to higher T c max by suppression of static stripe 
The Van Hove scenario |25|] finds no support in 
because it is the saddlepoint of the anti-bonding 
band which is at the LDA Fermi level in YBa2Cu307; 
the bonding band is about half-filled and enhances spin- 



fluctuations with q » (71", 7r) p6[ . The propensity to buck- 
ling is increased by pushing the conduction band to- 
wards the O a /bp z -level by lowering of e s but re- 
cent structural studies as well as Fig. 5, disprove 
that static buckling enhances T c max , although dynami- 
cal buckling might. The interlayer-pair-tunnelling mech- 
anism [^7| is ruled out by the fact that the additional 
factor cos ^k x cos ^k y attained by t (k) in bet materials 
suppresses the interlayer pair-tunnelling in Tl2Ba2CuOg 
compared with HgBa2Cu04, and yet, T cmax ^90 K in 
both cases. That the axial orbital is the channel for cou- 
pling the layer to its surroundings is supported [^8| by 
the observations that the k-dependence of the scattering 
in the normal state is w 2 -like || and that the c-axis trans- 
port is strongly suppressed by the opening of a pseudogap 
|p9| with similar k-dependence. The axial orbital is also 
the non-correlated vehicle for coupling between oxygens 
in the layer. Therefore it seems plausible that contraction 
of the axial orbital around the Cu02-layer, away from the 
less perfect doping and insulating layers, will strengthen 
the phase coherence and thus increase T cmax . Thermal 
excitation of nodal quasiparticles ||(J is, on the other 
hand, hardly the mechanism by which the superconduct- 
ing state is destroyed, because the axial orbital does not 
influence the band in the nodal direction. Finally, we 
note that the correlation between r and T c max does not 
extend to electron-doped cuprates, where the mechanism 
for superconductivity thus seems to be different. 

Discussions with H. Beck, I. Bozovic, and Z.-X. Shcn 
are gratefully acknowledged. 
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